Interaction-induced localization of anomalously-diffracting nonlinear waves 
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We study experimentally the interactions between normal solitons and tilted beams in glass waveg- 
uide arrays. We find that as a tilted beam, traversing away from a normally propagating soliton, 
coincides with the self-defocusing regime of the array, it can be refocused and routed back into 
any of the intermediate sites due to the interaction, as a function of the initial phase difference. 
Numerically, distinct parameter regimes exhibiting this behavior of the interaction are identified. 
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Solitary waves in general, and optical solitons in par- 
ticular, are the most important manifestations of the in- 
terplay between nonlinear self-focusing and the diffrac- 
tion/dispersion, which are inescapable characteristics of 
linear wave-packet propagation. Periodic structures, 
such as weakly-coupled waveguide arrays (WGAs) intro- 
duce two new important features to this problem. First, 
they support discrete solitons (DS), which are robust 
against collapse with increasing input power and prop- 
agation length Second, they allow control over the 
diffraction properties of tilted beams, giving rise to re- 
gions of normal, zero or anomalous diffraction for linear 
excitations p(, and accordingly, regions of self- focusing 
(SF) or self-defocusing (SDF) for nonlinear excitations 
y . These controllable characteristics not only make such 
structures favorable candidates for applications in opti- 
cal switching. They also make optical WGAs a model 
for a diverse array of classical and quantum systems 
that have nontrivial dispersion properties and support 
localized nonlinear modes, from Bose-Einstein conden- 
sates in periodic optical lattices 0, , to localized lat- 
tice vibrations pj and nonlinear excitations in arrays of 
mesoscopic- mechanical cantilevers [sj. 

Recently, there is a growing interest in interactions be- 
tween co-propagating beams in these WGAs 0, > m ~ 
eluding interactions between spatially-separate solitons 
[ill \vA H^ . as well as between co-propagating soli- 
tons and linear waves 

E M Gl M In homoge- 
neous Kerr media, mutually coherent in-phase and out- 
of-phase solitons attract or repel each other, respectively 
was found theoretically 9] and experimentally 
that the collision between a DS and a weak non- 
diffracting signal beam can lead to dragging of the blocker 
soliton towards the weak signal beam. In Ref. [12], 
Meier et. al. have experimentally studied the interac- 
tion between two narrow collinear DS in Kerr-nonlinear 
WGAs fabricated of AlGaAs. It was found that a pair 
of DSs, separated by 3 sites, are stable (i.e. noninter- 
acting) when the relative phase difference between them 
is around n, and unstable (i.e. attractive) when it is 
close to zero, in accordance with theoretical predictions 



FIG. 1: (Color online) Formation of a single normally injected 
soliton in the glass waveguide array, for a 40 /im wide input 
beam. The three panels correspond to different imaging con- 
ditions, and in each one the input power increases from bot- 
tom to top. (a) shows images of the output facet of the sam- 
ple with the imaging lens adjusted for the laser wavelength. 
At low powers the discrete-diffraction pattern is observed; At 
high powers a soliton forms, but is out of focus, (b) The soli- 
ton becomes visible at a different focal plane, corresponding 
to the Raman-shifted wavelength, but its image is saturated, 
(c) The same image as in (b), with an ND2 filter (X100 atten- 
uation), showing that the soliton is much stronger than the 
background. Setting (c) is used to study soliton interactions. 
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|10|,|16|- However, interactions involving tilted beams, 
at different diffraction regimes, have never been studied 
experimentally. 

In this paper, we report results for the interaction be- 
tween a normally-injected soliton (NS) and a nonlinear 
tilted beam (TB) that propagates away from the NS. We 
find that the characteristics of the interaction are consid- 
erably different in the SF and SDF regimes of TB prop- 
agation. Specifically, in the SDF regime we demonstrate 
for the first time that it is possible to restore self-focusing 
to the TB, and to redirect it into any of the intermediate 
sites, by varying the relative phase difference. 

We use a WGA fabricated of silica glass, as in Ref. 
[19], with the periodicity d = 12 fim and coupling con- 
stant C = 0.23 mm -1 . The sample is 13 mm long, cor- 
responding to 1.8 coupling lengths. As the Kerr coeffi- 
cient in glass is ~500 times smaller than in AlGaAs [TH . 
shorter pulses with higher peak powers are required to 
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FIG. 2: Experimental setup. OPA - optical parametric am- 
plifier, BS - beam splitter, ND - neutral density filter, CL - 
cylindrical lens, AL - aspherical lens, MO - microscope objec- 
tive. Inset: Lateral intensity profile of the input beams. The 
normal beam excites ~3-4 sites, while the tilted beam ex- 
cites ~4-5 sites. The excitation parameters, AO and K x , de- 
scribed in the text, are determined from applied translations 
and rotations, respectively, using the free-space wavenumber 
ko — 2tv/X . 



generate nonlinear effects. Importantly, the lowest-order 
nonlinear loss mechanism in glass is six-photon absorp- 
tion, which has an extremely low probability amplitude. 
DSs in glass WGAs, including TBs, are therefore sus- 
tained at high powers without breakup [19j |. allowing a 
broader range of powers and tilts to be explored. Glass is 
also characterized by strong stimulated Raman scatter- 
ing, which expresses itself as self frequency-shift of the 
soliton's spectrum to longer wavelengths [l8|. The for- 
mation of a soliton is therefore accompanied by a change 
of the focusing plane, due to the chromatic aberration 
of the output microscope objective lens. This sudden 
change of focus enables us to distinguish between strong 
solitary features and the linear background without ex- 
plicit spectral measurements, as demonstrated in Fig. 1. 

The experimental setup for the 2-beam experiment is 
sketched in Fig. 2. A train of 70-fs pulses at Ao=1.52 
/im, with < 20 MW peak power, was split into two iden- 
tical beams. The normally injected beam forming the NS 
was given an adjustable phase delay AO using a variable 
delay stage, and shaped using a cylindrical lens. The 
second beam passed through a rotating cylindrical lens 
that moved the beam on the aspherical coupling lens, 
and thus changed its angle of incidence a at the input 
facet of the sample. This angle determines the trans- 
verse wavenumber K x and thus the phase difference be- 
tween adjacent sites in the TB excitation. Both beams 
were coupled into the sample through the same aspheri- 
cal lens, which has a large numerical aperture (0.68 NA) 
and is optimized for the operating wavelength. A unique 
feature of this setup is the large clear aperture of the 
coupling lens, which in combination with its small focal 




FIG. 3: (Color online) (a),(b) Output profiles as a function 
of the phase difference, with both beams at normal incidence. 
The input powers are (a) 0.15 MW, and (b) 0.7 MW (Lin- 
ear and soliton regimes, respectively), (c) The output profiles 
of noninteracting 0.7 MW beams as a function of the TB 
tilt (data obtained with a very large time delay between the 
pulses). (d),(e) Interactions of 0.7 MW beams as a func- 
tion of the phase difference for the tilt angles of (d) 0.6° 
(K x = 0.17tt/cO, and (e) 3° (K x = OMn/d). The bottom 
parts in (a),(b),(d) and (e) show the output in the absence of 
interaction. 



length, results in a wide range of accessible incidence an- 
gles for the tilted beam, while maintaining paraxiality 
conditions and good coupling efficiency. The edge of the 
first Brillouin zone for this sample is attained at an in- 
cidence angle of a = arcsin(Ao/2<i) = 3.6°. In order to 
have sufficient resolution in K X1 that allows probing of 
the different diffraction regimes for the TB, the lateral 
width of the input TB had to be larger than ~30 jam. 
The combined input intensity profile for the two beams 
is shown in the inset of Fig. 2. The TB input profile did 
not vary under relevant tilts. 

Experimental results are shown in Fig. 3. Panels 
(a),(b) show the results obtained when the two beams 
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FIG. 4: (Color online) Results of numerical simulations, (a) Top: the amplitude distribution in the input, |<3> n (0)|. Bottom: 
output intensity distributions as a function of the tilt K x for noninteracting beams (i.e., superposition of individual solutions). 
(b),(c) Output intensity distribution as a function of the phase difference for: (b) K x — 0.247T and (c) K x — 0.97T. (d)-(i) 
Propagation maps of the intensity distribution as a function of z for < z < 13 mm and K x = 0.9tt: (d) NS without TB; (e) 
TB without NS; (f)-(i) both beams present, with the relative initial phase differences: (f) 0, (g) 7r/2, (h) 7r, and (i) 3tt/2. 



are at normal incidence. At low input powers (a), linear 
interference fringes are superimposed on the superposi- 
tion of the two overlapping discrete-diffraction patterns. 
It was found that the slope of the fringes exactly coin- 
cided with the phase slope imposed by the translation 
motor. At high input powers (b), the anticipated nonlin- 
ear instability around A6 — 2tt, which partially diverts 
the beams into intermediate sites, is observed |l2j . 

As tilts are imposed on the right DS (Fig. 3(c)), it is 
routed aside. As expected, the sideways shift of the beam 
increases in the SF regime, and decreases again in the 
SDF regime, where the beam becomes delocalized (The 
data in Fig. 3(c) was obtained with a very large time de- 
lay between the pulses, which inhibits the interaction). 
Panels (d) and (e) of Fig. 3 show examples of interacting 
NS and TB in the SF and SDF regimes of the TB, re- 
spectively. For small tilts that correspond to a TB in the 
SF regime (Fig. 3(d)), the observed attractive instabil- 
ity near 2tt phase multiples is less pronounced compared 
to the collinear excitation (Fig. 3(b)). This may be ex- 
pected due to the decreasing overlap of the two fields as 
the TB is routed aside. However, when the tilt angle is 
in the SDF regime of the TB (Fig. 3(e)), the instability 
is extended for phase differences far beyond a vicinity of 
2tt. Remarkably, the TB can be routed as a localized 
DS into any of the intermediate sites, with the appropri- 
ate phase setting, while the NS is partially annihilated. 
Evidently, the WGA favors interaction when one of the 
beams propagates in the SF regime and the other one in 
the SDF regime, as demonstrated by the extended range 
of phases exhibiting interaction dynamics. 

In order to gain better understanding of the interaction 
dynamics, we solved the one-dimensional discrete nonlin- 
ear Schrodinger equation in dimensionless form 0, [2(j , 
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h- 2 ($ n+1 + $ n _i - 2$ n ) + |$ n | 2 $ n = 0, (1) 



tween the characteristic nonlinear length Z^l and the 
diffraction length Zd = 1/C, is a single parameter 
characterizing the sample. The input is taken as two 
Gaussians centered at waveguide numbers n c \ = 51 and 

n c2 55, 
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on a grid of 101 sites; h 2 = ZnlC, the ratio be- 



where P* is the characteristic soliton power that de- 
fines Z NL Hj. We took P* = 0.5 MW. This value 
and the other experimental parameters (e.g. the diffrac- 
tion length, mode area, and Kerr coefficient) yield h~ 2 = 
4.9 x 10" 3 . The input widths, w\ = 2.5 and w 2 = 3.5, 
correspond to the experimental conditions. Finally, in 
these units the first band of linear waves corresponds to 

< K x < 7T. 

We have found that, to reproduce the propagation 
dynamics observed in the experiment, two conditions 
must hold regarding the excitation amplitudes. First, 
the TB amplitude, for the above input width, must be 
45 < Atb < 60. For higher amplitudes (with Ans = 0) 
the steering of the TB to nearby sites in the SF region 
is suppressed, and the soliton instead becomes trapped 
in one waveguide [HI EH- Simultaneously, it becomes 
unstable under small variations of the tilt. The second 
condition pertains to the NS. While the excitation am- 
plitude must be 40 < Ans < 90 for a NS to form by itself 
(i.e. with Atb = 0), we have found that it must be com- 
parable to the TB amplitude when the two excitations 
are sent together with the initial condition (2). Other- 
wise, the TB vanishes during the propagation, with only 
a small effect on the NS. In addition to the above con- 
straints on the excitation amplitudes, we have also found 
that the balance between the strength of diffraction and 
nonlinearity must be in favor of the nonlinearity, in the 
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form of a small h~ 2 parameter, to avoid breakup of the 
TB into filaments due to modulational instability [2^ . 
Indeed, all of the above conditions were fulfilled in our 
experiment. 

Results of the numerical simulations are shown in Fig. 
4, for A NS = 48 and A TB = 40. Figure 4(a) shows the 
amplitude profile of the input beam and the output inten- 
sity distribution for individual (non-interacting) beams, 
as a function of the TB tilt. Figures 4(b) and 4(c) show 
the output intensity distribution as a function of the ini- 
tial phase difference, for a TB that is SF and SDF, re- 
spectively. The calculated results are in good qualitative 
agreement with the experimental data displayed in Fig. 
3. In particular, Fig. 4(c) indeed shows routing of a relo- 
calized TB to nearby sites, which changes gradually over 
the entire range of the phase difference, < A6 < 2tt. 
The propagation maps as a function of z, panels (d)- 
(i), demonstrate how the NS is re-phasing the otherwise 
defocusing TB to form a localized soliton. The initial re- 
phasing generates a DS around one of the intermediate 
sites, which then periodically hops outwards. While the 
initial site in which the DS forms is a function of A#, the 
hopping period is always the same for a given K x , cor- 
responding to a constant propagation angle. The good 
agreement with the experimental data suggests that the 
spectral and temporal dynamics are not crucial for the 
understanding of the observed interactions. 

In conclusion, we have studied experimentally the in- 
teractions among two non-collinear discrete solitons in 
glass waveguide arrays. We have demonstrated for the 
first time that a tilted beam (TB) launched in a self- 
defocusing direction, while deviating away from its nor- 
mally propagating counterpart (NS), can be re- localized 
and routed into any intermediate site by the nonlinear 
interaction with the NS. The routing is controlled by the 
initial phase difference between the solitons. In simula- 
tions, this behavior was reproduced under several condi- 
tions, in which (i) the nonlinear length is considerably 
smaller than the diffraction length; (ii) the NS and TB 
have comparable input amplitudes; and (hi) the TB in- 
put amplitude is moderate, to enable soliton steering and 
to avoid its trapping and instability. When two on-site 
excitations exist simultaneously, <£ n = A + B, the non- 
linear Kerr term in the coupled equations (1) acts as a 
non-trivial interference term |<l> n | 2 <l> n = (\A\ 2 + \B\ 2 + 
2\A\\B\sinA6)(A + £?), and a new propagation mode 
that is in balance with the weak inter-site coupling is 
obtained. Under the above conditions for the initial col- 
lective excitation, this mode favors localization of the 
initially SDF TB, which can indeed be controlled with 
the initial phase difference as a single parameter. The 
result is a unique recapturing of a defocusing beam into 
a soliton, which should have analogs in many other phys- 
ical systems that have nontrivial dispersion properties 
and support localized nonlinear modes. In addition, the 
interaction that we observe has potential applications in 



all-optical switching devices. 
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